Analytical and Chemometric Characterization of Sweet Pedro Ximénez Sherry Wine during Its Aging in a Criaderas y Solera System

Pedro Ximénez is a naturally sweet sherry wine produced in southern Spain from raisined Pedro Ximénez grape must and aged using a traditional Criaderas y Solera system. Complete analytical characterization has been useful in determining which parameters are the most influential in the aging of this wine. The organic acids, volatile compounds (higher alcohols, esters, aldehydes, and acetals), and phenolic compounds of this wine evolve during its aging, mainly through physico-chemical reactions and the contributions of wood compounds. During their aging, Pedro Ximénez sherry wines develop their organoleptic profiles, as tasting sessions have confirmed. A strong correlation between the aging of a wine and the parameters analyzed has also been corroborated through an MLR analysis. This allowed for the development of a model that, by using just 8 of the variables considered in the study, led to the determination of wine samples’ ages at over 97% confidence. This constitutes a rather useful tool for wineries to control Pedro Ximénez sherry wine aging processes.


Introduction
Amontillado, Fino, Oloroso, Palo Cortado, and certain natural sweet wines, such as Moscatel or Pedro Ximénez, are produced in the Jerez-Xérès-Sherry protected designation of origin (PDO). The dry fortified wines are made from Palomino grapes, while the natural sweet wines are made from their respective varieties of Moscatel or Pedro Ximénez grapes. Aging is carried out according to a traditional dynamic aging system from the sherry area known as Criaderas y Solera. All of this, along with the specific climate conditions in Marco de Jerez, make these wines a highly appreciated product by oenologists worldwide [1,2].
Pedro Ximénez sherry wine (PXSW) is a liqueur wine with a protected designation of origin (Jerez-Xérès-Sherry). Specifically, it is a natural sweet wine obtained by fortifying alcohol, including alcohol, distillate, or brandy, with overripe or raisined grape musts of the Pedro Ximenez variety, in addition to an oxidative aging process for at least two years [3]. The freshly harvested overripe grapes at over 13.5 • Baumè are directly exposed to sunlight, on esparto grass or mesh nets in specifically appointed areas known as "paseras", for a period of between 5 and 15 days, depending on weather conditions. Thus, partial dehydration and darkening of the berries take place, which result in raisins with a reducing sugar content over 450 g/L. The raisins are then pressed to obtain a very sweet must, which is then fortified with wine alcohol up to between 15 and 17% ABV (alcohol by volume) in order to stop the fermentation process [4][5][6]. Young Pedro Ximénez fortified wine (YPXFW) is aged in a winery through oxidative aging under a Criaderas y Solera system. It then Foods 2023, 12,1911 2 of 20 acquires a more or less intense ebony color with iodized iridescences and a very dense visual appearance (it is probably the sweetest wine in the world) [4]. It is very complex in aromas, where notes of dried fruits, such as raisins, figs, or dates, stand out, together with others of honey, syrup, and candied fruits, with long-aged ones beginning to develop hints of toast (coffee, bitter chocolate, and cocoa) or liquorice. It possesses a natural acidity that comes from the raisined grapes, which balances its sweetness and makes of it a fresh, harmonious, velvety, and unctuous wine on the palate, with a very pleasant finish [4][5][6].
PXSW ages under a dynamic system of Criaderas y Solera, characteristic of sherry wines. This system consists of oak casks of 500-600 L of capacity, which are arranged at several scales or levels, according to the aging time of the wine that they contain [7,8].
The quality and botanical origin of the oak casks that contain aging PXSW, as well as the treatment and previous use that they have had, are of great relevance with regard to the aging process (these casks must have been previously seasoned before being used as part of a Criaderas y Solera system [9]). Casks, themselves, are not inert elements, but rather contribute with numerous compounds from their own compositions to aging wines and facilitate the occurrence of specific chemical reactions in them that have an impact on wines' organoleptic properties [10][11][12]. Wood mainly contributes with phenolic compounds [13], especially ellagitannins and phenolic aldehydes derived from cinnamic acid and benzoic acid, such as coniferaldehyde, sinapaldehyde, vanillin, and syringaldehyde [13,14]. Additionally, furfurals and their derivatives, such as 5-methylfurfural and 5-hydroxymethylfurfural, are obtained from the wood [14]. Other substances that are transferred by the wood of casks into the wine during its aging are pentoses, polysaccharides, lactones, fatty acids, inorganic substances, and alcohols, which confer complexity and quality to the products [14].
In this work, complete physicochemical characterization of PXSWs aged through a Criaderas y Solera system as well as a sensory study of these wines with an average age of 40 years in the Solera have been carried out. A chemometric study has also been conducted in order to determine which variables might present any relevant correlations between the parameters considered and the aging times of the wines that would lead to the development of multiple linear regression (MLR) models that would allow for the estimation of the aging time of the PXSWs simply by identifying the most influential parameters involved in their aging processes. This model should represent a useful tool for wineries to regularly control the quality of their wines.

Samples
Bodegas Fundador, S.L.U. (Jerez de la Frontera, Cadiz, Spain) provided the PXSWs, the casks, the wine distillates for the fortification used, and the winery facilities where this study was conducted. The casks used in this study have been selected from the industrial Criaderas y Solera system at Bodegas Fundador, S.L.U, and are old enough for us to consider that the system of continuous mixtures on which the Criaderas y Solera systems are based have reached equilibrium on all scales. All of them were made of American oak (Quercus alba), with medium wood grain, medium-toasted, and 500 L of capacity. Figure 1 displays the process diagram being studied. All wine replenishments for the different aging scales were carried out using the PXSW from the next younger aging scale previously fortified at 18% ABV, except for the 5th Criadera, which was replenished with YPXFW at 17.5% ABV. This study lasted four years, and each wine category was sampled during the mo of May of each year, as described below: -Young Pedro Ximénez fortified wine (YPXFW) (vintage): Each year a combined sa ple was taken from 5 industrial 20,000 L tanks containing the young fortified w from that year's harvest. One liter was taken from each of the tanks and blended obtain a final five-liter sample. Before proceeding with the sampling, the wine fr the tanks was tasted to ensure that it presented no organoleptic flaws. -Pedro Ximénez sherry wine (PXSW) aged in a Criaderas y Solera system: The sys consisted of one Solera and five Criaderas that comprised 25 × 500 L casks per s and involved the performing of regular removals, known as sacas, and replen ments, known as rocí os. The 5th Criadera was an average of 2 years old (YO); the Criadera, 6 YO; the 3rd Criadera, 12 YO (PX12YO); the 2nd Criadera, 20 YO; the Criadera, 30 YO; and the Solera, 40 YO (PX+30YO). Each year, over the whole 4-y research period, a combined 5 L sample was gathered by extracting 200 mL of w from each of the 25 casks that made up each of the scales. Before proceeding with sampling, the wines from each of the casks were tasted in order to verify that organoleptic flaws were present.
A combined sample from each scale (6 different scales + YPXSW) taken over th years of the study, i.e., 4 samples from each wine age, was physicochemically charac ized. All the samples were analyzed in triplicate. Only the samples taken in the fou year were used for the sensory evaluation.

Reagents
Ultrapure water (EMD Millipore, Bedford, MA, USA); UHPLC-grade acetone (V International, Radnor, PA, USA); and 0.1 M sulfuric acid (Sigma-Aldrich, Saint Louis, M USA) were used for the preparation of the eluent used for the analysis of organic acid Folin-Ciocalteau reagent and anhydrous sodium carbonate from Merck (Darmst Germany), as well as ultrapure water (EMD Millipore, Bedford, MA, USA), were used determine the Folin-Ciocalteau Index.  This study lasted four years, and each wine category was sampled during the month of May of each year, as described below: -Young Pedro Ximénez fortified wine (YPXFW) (vintage): Each year a combined sample was taken from 5 industrial 20,000 L tanks containing the young fortified wine from that year's harvest. One liter was taken from each of the tanks and blended to obtain a final five-liter sample. Before proceeding with the sampling, the wine from the tanks was tasted to ensure that it presented no organoleptic flaws. -Pedro Ximénez sherry wine (PXSW) aged in a Criaderas y Solera system: The system consisted of one Solera and five Criaderas that comprised 25 × 500 L casks per scale and involved the performing of regular removals, known as sacas, and replenishments, known as rocíos. The 5th Criadera was an average of 2 years old (YO); the 4th Criadera, 6 YO; the 3rd Criadera, 12 YO (PX12YO); the 2nd Criadera, 20 YO; the 1st Criadera, 30 YO; and the Solera, 40 YO (PX+30YO). Each year, over the whole 4-year research period, a combined 5 L sample was gathered by extracting 200 mL of wine from each of the 25 casks that made up each of the scales. Before proceeding with the sampling, the wines from each of the casks were tasted in order to verify that no organoleptic flaws were present.
A combined sample from each scale (6 different scales + YPXSW) taken over the 4 years of the study, i.e., 4 samples from each wine age, was physicochemically characterized. All the samples were analyzed in triplicate. Only the samples taken in the fourth year were used for the sensory evaluation.

Reagents
Ultrapure water (EMD Millipore, Bedford, MA, USA); UHPLC-grade acetone (VWR International, Radnor, PA, USA); and 0.1 M sulfuric acid (Sigma-Aldrich, Saint Louis, MO, USA) were used for the preparation of the eluent used for the analysis of organic acids.
Folin-Ciocalteau reagent and anhydrous sodium carbonate from Merck (Darmstadt, Germany), as well as ultrapure water (EMD Millipore, Bedford, MA, USA), were used to determine the Folin-Ciocalteau Index.
For the analysis of phenolic and furfural compounds, HPLC-grade acetonitrile, supplied by Panreac (Barcelona, Spain); acetic acid, supplied by Merck (Darmstadt, Germany); and ultrapure water, supplied by EMD Millipore (Bedford, MA, USA) were used to prepare the UHPLC phases.
Sigma-Aldrich (Saint Louis, MO, USA) supplied all of the reagents for the analysis of the oenological control parameters in addition to the calibration standards.

Parameters of Oenological Control
The official methodology described by the International Organisation of Vine and Wine (OIV) was used to determine the oenological control parameters: The density of wines (g/L) was directly measured by means of a DMA-5000 digital density meter (Anton Paar, Ashland, OR, USA). The alcoholic strengths (% ABV) were determined by distilling the wine and then measuring the distillate density [15] by means of a DMA-5000 digital density meter (Anton Paar, Ashland, OR, USA). The pH of wines was determined with a BASIC 20 pH meter (Crison Instruments SA, Barcelona, Spain). The total acidity of wines (g tartaric acid/L) was determined via potentiometric titrations at pH 7 [16]. The volatile acidity of wines (g acetic acid/L) was analyzed through the use of an AA3 HR Autoanalyzer segmented flow analyzer (Seal Analytical, Norderstedt Stadt, Germany) [17]. The glycerol content (g/L) was determined via an enzymatic method [18]. The total sulfur (mg/L) was analyzed through the use of the Ripper method [19]. The sulfate concentration (in g of K 2 SO 4 /L) was calculated via the gravimetry of BaSO 4 precipitation [20]. The calcium, potassium, copper, and iron contents (mg/L) were determined via atomic absorption spectroscopy, using a PinAAcle 900F system equipped with the software application WinLab32 AA, both by Perkin Elmer (Boston, MA, USA). The total dry extract of wines (g/L) was determined via gravimetry, following the official standard procedure [21]. The reducing substances (g/L) were also analyzed through following the official standard procedure [22]. The sugar-free extract of wines (g/L) was determined according to the following formula [21]: Sugar − f ree extract g L = Total dry extract g L − Reducing substances g L

Organic Acids
A 930 Compact IC Flex, from Metrohm (Madrid, Spain), with a Metrosep Organic Acids column of 9 µm particle size and 250 × 7.8 mm (i.d.) [8], was employed to analyze citric, lactic, malic, succinic, and tartaric acids (mg/L). The software application used for data acquisition and processing was MagicNet 3.3 (Metrohm, Madrid, Spain). The compounds were identified via a comparison of their retention time against the standard used.

Higher Alcohols, Aldehydes, Methanol, and Ethyl Esters
An Agilent 7890B Gas Chromatograph (Agilent Technologies, Santa Clara, CA, USA), coupled with a flame ionization detector [8,23], was used to determine methanol, acetaldehyde, acetaldehyde-diethylacetal, esters (diethyl succinate, ethyl lactate, ethyl acetate, ethyl hexanoate, ethyl octanoate, ethyl decanoate, ethyl dodecanoate, and ethyl tetradecanoate), and higher alcohols (isobutanol, n-propanol, 2-methyl-1-butanol, 3-methyl-1-butanol, nhexanol, and 2-phenylethanol). Standards were prepared in an ethanol/ultrapure water solution at the same alcoholic strength as the wines for the identification, based on their retention times, and quantification of the samples. For the ISTD (internal standard) for aldehydes and higher alcohols, 2-pentanol was used, and ethyl undecanoate was used for major esters. The samples were injected after distilling them at the same alcoholic strength.
The results have been expressed as mg/L.

Folin-Ciocalteau Index
The Folin-Ciocalteau index (FCI) was used to determine the total content in phenolic compounds [24]. A Lambda 25 spectrophotometer from Perkin Elmer (Boston, MA, USA) was used to measure absorbance at 750 nm on glass cuvettes with a light path length of 10 mm. Previously, a calibration curve was performed with gallic acid in a concentration range between 0 mg/L and 750 mg/L. Gallic acid standards were prepared in a solution at the same alcoholic strength as the wines, containing 425 g/L of a mixture of glucose and fructose in a 0.93:1.00 ratio, similar to that of wines. The results were expressed as mg gallic acid equivalent (GAE)/L.

Phenolic Compounds and Furfurals
Fifteen phenolic compounds (caffeic acid, p-coumaric acid, ferulic acid, gallic acid, p-hydroxybenzoic acid, protocatechuic acid, syringic acid, vanillic acid, trans-caftaric acid, cis-p-coutaric acid, trans-p-coutaric acid, fertaric acid, p-hydroxybenzaldehyde, syringaldehyde, and vanillin) and two furfurals (5-hydroxymethylfurfural and furfural) were determined via UHPLC [23,25]. A Waters Acquity UPLC fitted with an Acquity UPLC C18 BEH column, 100 × 2.1 mm (i.d.) and 1.7 µm particle size, as well as a PDA detector (Waters Corporation, Milford, MA, USA) has been used. Nylon membranes with a pore size of 0.22 µm were used to filter the samples and standards. Standards were prepared in an ethanol/ultrapure water solution at the same alcoholic strength as the wines. The analyzed compounds were identified via a comparison of the retention time and the UV-Vis spectrum of the sample, as well as with the standard used. The results were expressed as mg/L [23,25].

Determination of the Brown Color of the Wines
The measure of the absorbance at 470 nm was used to determine the brown color of the wines [26]. All of the measurements were carried out through the use of a Perkin Elmer spectrophotometer-a Lambda 25 spectrophotometer (Perkin Elmer, Boston, MA, USA). All of the results were expressed as absorbance units. Prior dilutions were made via the use of a hydroalcoholic mixture at 17% ABV (deionized water and 96% ABV neutral wine alcohol) when necessary. The measurements were performed in triplicate.

Tasting Sessions
A properly isolated room, in order to allow the tasters to focus on their job, was used for the tasting sessions. Individual tastings sessions were carried out at an ambient temperature of 22 • C [27]. The tasting panel was made up of 7 tasters from Bodegas Fundador's staff, and all of them had over 10 years' experience in sherry wine tasting.
Two PXSWs of different average ages were evaluated in the tasting sessions: 3rd Criadera (12 years), named PX12YO, and Solera (40 years), named PX+30YO, given that the specifications for the Jerez-Xérès-Sherry designation of origin (DO) only allow the certification of 12-year-old wines and over 30-year-old wines (VORS-very old rare sherry), as well as YPXFWs.
For each wine, 50 mL was served in black wine-tasting glasses [28], and these were capped, for at least ten minutes, with a glass lid in order to stabilize their headspace prior to tasting. Nose and mouth perception data were recorded.
The definitions associated with the descriptors, selected according to the criteria of the DO tasting panel, are specified in Table 1. The olfactory-gustatory standard that corresponds to the maximum intensity of most of them but dried fruit (value: 9) on the 9-point interval numeric scale used [29] was a 50-year-old PXSW (special casks). A young vintage wine obtained from the must of sunlight-unexposed overripe Pedro Ximénez grapes, to which distilled alcohol at 95.4% ABV was added to stop its fermentation (17% ABV and 217 g sugars/L), was used as the minimum intensity standard (value: 1) except for the descriptor dried fruits, where it showed the highest intensity [29]. Table 1. Olfactory descriptors and patterns used during the work with the tasters.

Descriptor Definition
Odor Aromatic intensity Intensity of all of the positive aromatic notes of the wine.
Dried fruits Aromas of ripe and partially dehydrated fruits due to direct exposure to sunlight (raisins, figs, dates, prunes, etc.).

Oak
Aromas of wood, with hints of vanilla and spices.
Flavor Sweetness Very sweet sensation but not sickening.

Equilibrium
Good integration of the alcohol, acidity, and sugar, without astringency, but with the aromatic reminder of the oak, as appropriate for an oak-aged wine.

Persistence
Time evaluation of the olfactory-gustatory notes remaining after the final sip.

Statistical Analysis
The Statgraphics 19 software package (Statgraphics Technologies, Inc., The Plains, VA, USA) was employed for the multiple linear regression analysis, ANOVA, and Fisher's least significant difference test. Microsoft Excel 2016 (Microsoft Corp., Redmond, WA, USA) was employed to process the rest of the statistical parameters.
Regarding the sensory data, the statistical analysis included an ANOVA, as well as a factorial analysis. The software application Statistica 8.0 (StatSoft Inc., Tulsa, OK, USA) was used for both of these analyses. Microsoft Excel 2016 (Microsoft Corp., Redmond, WA, USA) was used to generate a spider web chart.

Alcoholic Strength
The alcoholic content of the PXSWs studied (Table 2) revealed a decreasing trend over their aging, where significant differences (ANOVA) were found corresponding to the different aging scales, since the high concentration of sugars in this wine prevents the staves that form the casks from tightly binding together and preventing ethanol losses via evaporation. This is due to the high density of this type of sweet wine as well as to the hygroscopicity of the sugars, which prevent the wine from sufficiently penetrating the wood pores and making them swell, as occurs in the case of the oxidative aging of dry fortified sherry wines (Amontillado, Oloroso, or Palo Cortado). In fact, in these cases, the trend is quite the opposite, with the alcohol content increasing as the wine ages [7,8,30,31]. These losses in ethanol and other volatiles, together with other minor losses attributable to the transpiration of the water molecules through the pores of the wood that form the casks containing the wine, is a phenomenon traditionally known as "merma" [8,32].
The average annual value corresponding to merma over the 4 years of the study was approximately 2.5% in the Criaderas y Solera system analyzed, i.e., a yearly reduction in the wine volume in the casks by 2.5%, which means an increment of the concentrations of the compounds present in the casks [7,8]. Mean values ± standard deviations (n = 4) are shown. ANOVA: for the same wine, different letters (in a line) indicate significant differences (p < 0.05). Cra: Criadera; ac.: acid. (a) Density (g/L); alcoholic strength (% ABV); total acidity (g TH 2 /L); volatile acidity (g AcH/L); and glycerol (g/L); (b) organic acids (mg/L); sulfates (g K 2 SO 4 /L); calcium; potassium and total SO 2 (mg/L); total dry ext. (total dry extracts, g/L), red. subst.

Acidic Parameters
The ANOVA results do not reveal a clear trend in the wines' pH values; significant differences were found between the total acidity corresponding to the different aging scales, as well as between their volatile acidity levels ( Table 2).
A decrease in pH with aging (4.50-4.06) was observed in the PXSW Criaderas y Solera system ( Table 2). Plastering is not used for the vinification of Pedro Ximénez, nor are corrections of the fixed acidity with tartaric acid at harvesting employed, which result in a higher pH and a lower total acidity in relation to other sherry wines where these practices are applied [7,8,30,33,34]. The stabilization of pH at around 4.0 in the PXSW from the 1st Criadera and from the Solera must be a consequence of the increased concentration caused by the merma of the acids from the raisined grapes and to the precipitation of insoluble salts that seek a chemical equilibrium by concentrating their ions (mainly potassium bitartrate and calcium tartrate).
The total acidity in the PXSW Criaderas y Solera system increases with the aging time (Table 2), from 2.95 g tartaric acid/L in the YPXFW to 5.25 g tartaric acid/L in the Solera PXSW. This is explained by several factors, such as the following: the increases in acetic acid content via the oxidation of the ethanol in the medium [23]; the concentration effect resulting from the merma of organic and inorganic acids in the wine; and the contributions of acidic substances from the wood (mainly acetic acid [14,35] and fatty acids [14]). This is the case despite the decrease in the concentration of tartaric acid that takes place in the wine during its aging as a consequence of the precipitation of its potassium and calcium salts.
Likewise, an increase in volatile acidity during aging is observed in these wines, from 0.52 g acetic acid/L in the YPXFW to 1.06 g acetic acid/L in the Solera PXSW (Table 2). This is due to the following factors: the oxidation of the alcohol in the medium [23]; the contributions of acetic acid from the wood [14]; and the effect of the merma. In addition, the volatile acidity of solely the Solera wine shows that the product has been aged correctly from a microbiological point of view (without bacterial infections); this value is much lower Foods 2023, 12,1911 8 of 20 than 2.10 g acetic acid/L, which is the maximum value according to the specifications set out for the Jerez-Xérès-Sherry DO corresponding to this type of old wine [36].

Glycerol
Significant differences (ANOVA) have been found in the glycerol content because of aging for all of the wines studied. The evolution of glycerol concentration was upward with a longer aging time (Table 2), and reached values close to 10 g/L in the Solera wine. The concentration of this compound is caused by the evaporation and/or transpiration associated with the merma process.
In wine, glycerol is mainly a secondary product of alcoholic fermentation [37], although it is also explained by the hydric stress suffered by grape cells during their raisining. This process triggers an intracellular anaerobic metabolic process that leads to the synthesis of glycerol [38], which acts as an osmolyte. In the young PXSW, despite the discontinuation of the fermentation process and because of the two effects mentioned above, some wines with contents around 3 g/L glycerol have been reported.

Sulfates
Although the differences between the first and second scales with respect to sulfate content were not significant, as the aging time grew longer statistically significant differences (ANOVA) between the scales could be observed. The presence of sulfates is exclusively attributable to the grapes, but an increment was observed according to aging time (Table 2), mainly due to the effect of the merma. Nevertheless, although a greater rise has been expected, this did not take place, which could be due to precipitations, mainly of calcium sulfate [39].

Potassium and Calcium
Significant differences (ANOVA) have been found in most of the wines studied between potassium and calcium contents as a function of aging [40,41]. In the younger scales of the aging system, a decrease in potassium and calcium content was observed (Table 2) with respect to the initial wine, while in the older wines their content increased due to the merma. This reduction in the younger wines was explained by the precipitation of inorganic and organic salts that takes place during their aging. In the wines' hydroalcoholic medium, potassium bitartrate is supersaturated, which makes it precipitate during the aging in an attempt to achieve an equilibrium between potassium and tartaric acid. These precipitations are favored by the temperature drops that occur in cellars during the winter. In the case of calcium, it is the calcium salts of tartrate and sulfate that precipitate. In certain wines with a long aging period [40,41], oxalate also precipitates because of the higher-than-usual concentration that results from the merma.

Total Sulfur Dioxide
Significant differences (ANOVA) could be observed between the values corresponding to total sulfur dioxide in the different aging scales, since no new sulfur was added after the crushing and sulfur oxidizes into sulfate. The values registered for the YPXFW were a consequence of the additions that take place during the crushing of the grapes (90-120 mg/L). On the other hand, during aging, a drop in the total sulfur content of the wine, which reached almost 0 mg/L in the Solera, was observed ( Table 2).
3.1.7. Density, Total Dry Extract, Reducing Substances, and Sugar-Free Extract Significant variations (ANOVA) were observed between the levels of these parameters (total dry extract, reducing substances, and sugar-free extract) and the different aging scales.
The density value of PXSWs is much higher than that of other types of sherry wines [7,8,30], due to their high sugar contents. This parameter also evolved slightly upward during the aging studied, and significant discrepancies have been observed between the density values corresponding to the wines from the different aging scales.
During the aging, an increment in the content of soluble solids as well as an increasing density of the wines due to the merma were observed. The poor air tightness of the casks allows for alcohol losses that favor this upward tendency in soluble solids and density.
Total dry extract, reducing substances, and sugar-free extract are parameters that increased during the aging process of the PXSWs studied.
Thus, the total dry extract increased with aging due to the merma phenomenon, which affects all of the compounds involved in its quantification, i.e., acids, glycerol, polyphenols, and reducing substances, such as the hexoses, pentoses, and polysaccharides extracted from the wood or the glucose and fructose originating from the Pedro Ximénez raisins. PXSWs have much higher total dry extract values than the dry sherry wines [7,8,30] because of the greater amount of sugars in them.
With regard to the oxidative aging, the concentration of reducing substances as a consequence of the merma and also because of the contributions from the wood's hemicellulose [13,14] increased with aging time. This parameter is mainly associated with the quantification of sugars and polysaccharides, whose content is greater in PXSWs than in other sherry wines because of the presence of the glucose and fructose originating from the raisins [7,8,30].
Sugar-free extract is a measure of all non-volatile substances at 110 • C. The origin of these substances can be varied, including the following: the raisins themselves; the different fermentation reactions that take place in the wine before it is fortified to stop the alcoholic fermentation process; the wood of the cask where the wine is aging; or the merma effect, among others. This analytical parameter rose with oxidative aging from 20.83 g/L in the YPXFW up to 31.03 g/L in the Solera wine.

Organic Acids
Significant differences (ANOVA) were observed between the concentrations of each organic acid depending on the scale in the system studied. The concentrations of lactic, malic, citric, succinic, and tartaric acids in PXSWs depend mainly on their content levels in the grapes themselves, i.e., on the anaerobic metabolic reactions that take place in the cells of the grape berries during their raisining process [42,43], as well as on the greater or lesser activity that yeasts and bacteria may have had between the crushing of the raisins and the fortification of the must, since the presence of these microorganisms in their metabolic pathways may decrease or increase in the wines [37]. Tartaric acid presented a tendency to decrease during the aging of the wine because of the insolubilization of its calcium and potassium salts [8,39]. The other acids increased their presence (Table 2) as a result of the merma.
The main sources of citric acid in YPXFW are raisins and stems [44]. The different acids in grape stems are usually salified. Citric acid, therefore, reaches the young wine through the pressing of the raisin bunches, where the presence of the stems enhances the extraction of the must from inside the raisins. It can be seen in Table 2 that, during the aging, there was an increment of citric acid from 210 to 385 mg/L. A similar behavior is observed for malic acid, and its origin is also the raisined grapes [42,43]. In the absence of either any significant metabolic activity by bacteria or yeasts or precipitations in which this compound may intervene, it becomes the organic acid at the highest concentration in the oldest wines, reaching almost 2.3 g/L in the Solera ( Table 2). The concentration of this acid increases during the aging process as a result of the merma.
The metabolic activity of yeasts during must fermentation produced succinic and lactic acids [37]; this is why their respective contents are not very high in PXSWs, where the fermentation process is halted just after it has started. A growing concentration, owing to the merma, is also observed as the wine is aged (Table 2); however, their levels at 514 mg succinic acid/L and at 653 mg lactic acid/L in the Solera wine are lower than those found in other Sherry wines, where fermentation is completed before they are aged under a similar oxidative aging system [7,8,37].

Aldehydes, Methanol, and Higher Alcohols
Significant dissimilarities (ANOVA) between the different aging scales were observed with respect to these volatile compounds, except in those cases where no notable evolution was registered with the aging, such as n-hexanol.
Aldehyde content (acetaldehyde and its diethyl acetal) in YPXFW is not very high, and it depends on the amounts produced in the initial fermentation of the must, as well as on the amounts contributed by the 95.4% ABV wine distillate used for its fortification. Aldehyde content presents an upward trend during the aging (Table 3) as a result of alcohol oxidation; however, their content does not reach the expected levels in PXSWs, which may be explained by the poor air tightness of the casks in comparison to that of casks used to age dry sherry wines. This lack of tightness favors the loss of these highly volatile compounds through evaporation. The methanol content present in YPXFW (Table 3) is in accordance with the degradation level of the grape cell walls that takes place during the raisining process. Methanol is thus released due to the demethoxylation of the pectins catalyzed by pectin esterases [41,45]; additionally, the high pressure used for the crushing of the raisins causes the extraction of precursors of this compound from the solid parts of the grape bunches. The volatile composition of the wine distillates used for the fortification of the must also plays a relevant role. The trend that has been observed during the aging process is not clearly defined, such that the content levels remain more or less steady, probably because of two mutually compensating effects: evaporation because of the poor air tightness of the casks and the merma.
All of the higher alcohols studied evolved upward during aging (Table 3), mainly due to the concentration of these compounds with the merma effect. Alcoholic fermentation produces higher alcohols via different metabolic pathways that correspond to the decarboxylation/deamination of amino acids [37,39], as well as from the smaller contribution from the 95.4% ABV wine distillate used for the fortification of YPXFWs. Since alcoholic fermentation is minimal in YPXFWs, the concentrations of n-propanol, isobutanol, 2-methyl-1-butanol-1, and 3-methyl-1-butanol are lower than in other oxidatively aged sherry wines [7,8,31] where the must is fully fermented.
The n-hexanol content was low in all of the wines (Table 3), ranging from 0.22 to 0.32 mg/L, and remained stable throughout their oxidative aging.
During alcoholic fermentation, 2-phenylethanol is a compound produced by the yeasts from 2-phenylalanine [45]. Furthermore, its content increases during the aging of the wine due to either the hydrolysis of 2-phenylethyl acetate or to the oxidation of ethanol and its subsequent condensation with the phenols from the wood [46], as well as to the merma. Its contents, shown in Table 3, are lower than those in other fully fermented and oxidatively aged sherry wines [7,8].

Major Ethyl Esters
Significant differences (ANOVA) in major ethyl esters were observed depending on the system scales. Among the esters studied (Table 3), ethyl acetate stands out as one of the compounds that most evolves and in an increasing manner with the aging of the wines analyzed. In the case of PXSWs, acetic acid is produced via the oxidation of ethanol in the medium [23], and there is also a contribution of acetic acid from the casks' wood [14]. Ethyl acetate is formed via the esterification of acetic acid with the ethanol in the medium, and it also increases its concentration because of the merma. Its maximum level is reached in the Solera wine at 129.0 mg/L, which is much lower than that registered for other oxidative aging Sherry wines after the same time in the aging casks [14]. This is partly due to the lower alcohol content of PXSWs compared to similarly aged dry sherry wines, which affects the equilibrium required to form the corresponding ester, but is also due to the aforementioned poorer tightness of the casks containing Pedro Ximénez wines. Table 3 also shows the results corresponding to other ethyl esters. It should be noted that the esters studied belong to two families: The first one, the fatty acid ethyl esters (ethyl hexanoate, ethyl octanoate, ethyl decanoate, ethyl dodecanoate, and ethyl tetradecanoate), are partially produced by the esterification with ethanol of the fatty acids produced by yeasts during the glyceropyruvic fermentation of the must sugars and represent one of the most important groups of aromatic compounds that play a significant role regarding the organoleptic characteristics of wine [47]. These esters have also been detected in significant amounts in the musts from grapes that have undergone intracellular anaerobic metabolism (a process that is favored by grape raisining) [38]. The concentrations of these esters hardly reached 1 mg/L in the system analyzed, with ethyl octanoate and ethyl decanoate being the major ones. Moreover, only ethyl octanoate shows a growing trend with aging, which explains why no significant variations between the different aging scales could be observed in most cases. The concentration caused by the merma did not seem to affect the concentration of these esters in the old wines, since they are highly insoluble in media below 45% ABV. Regarding the contents of fatty acid ethyl esters, a similar behavior to that observed in other types of sherry wines was registered [7,8].
The second of the ester families are the ethyl esters of organic acids (ethyl lactate and diethyl succinate), which are generated via the chemical esterification of ethanol with their organic acids that come from the activity of the yeasts during the initial alcoholic fermentation that takes place in the must when making YPXFWs. Table 3 shows the concentrations measured in the Criaderas y Solera system and in the young wine. Thus, an increasing trend associated with aging time was observed (reaching 51.37 mg/L ethyl lactate and 11.17 mg/L diethyl succinate in the Solera).

Folin-Ciocalteau Index, Phenolic Composition, and Lignin-Derived Compounds
For the FCI and the phenols analyzed, significant differences (ANOVA) were found between the different aging scales. In the system studied, a constant increase in the FCI was observed with aging (Table 4). This was due to the transfer of compounds from the wood, the effect of merma, and the different chemical reactions that take place during the aging process, as well as those that involve the phenolic substances contributed by the grapes [14,48]. Total phenol contents are much higher in PXSWs than in the other sherry wines [7,8] because of the significant phenolic and furanic aldehyde concentrations resulting from the raisining of the grapes exposed to direct sunlight [4][5][6]49]. The PXSW Solera had an FCI of 2.871 mg GAE/L.
It was also observed that aging PXSWs (Table 4) contain certain lignin-derived compounds from oak wood, but their concentrations are lower than those found in dry sherry wines [7,8]. This may be due to the poorer penetration that these sweet wines achieved into the wood pores because of the high sugar contents and high viscosity in comparison to dry sherry wines. These, with a higher fluidity, as well as alcohol content, penetrate the wood more easily and achieve greater extractions of wood compounds.
Regarding the phenols that are specific to grapes (Table 4), it was observed that, as the aging time increases, the wines present a lower content of tartaric acid derivatives. This was possibly due to hydrolytic processes, since their free phenolic structure increased at the time of hydrolysis (p-coumaric acid, caffeic acid, and ferulic acid) [44], although a decrease in the contents of caffeic acid and ferulic acid was also observed later in lengthily aged wines, with the exception of p-coumaric acid. Furthermore, when analyzing the phenols derived from tartaric acid in the starting wines, it was noted that the sunlight exposure process undergone by Pedro Ximénez grapes caused hydrolysis and a decrease in their concentration in the YPXFWs, compared to those determined for wines of the Palomino variety [7]. Other phenols, such as gallic acid-contributed by the raisined grapes themselves as well as by the degradation of gallotannins yielded by the wood-or syringic acid, which comes in small quantities from the raisined grapes as well as, in greater amounts, from the wood lignin or the oxidation of its aldehyde (syringaldehyde) [13], increased their concentration as a consequence of the merma effect [50].
There are higher contents of p-hydroxybenzoic acid and protocatechuic acid in YPXFWs, because of the raisining effects [4,5,49], compared to dry sherry wines [7,8]. During aging, these contents remain at a high level and, in addition, they may increase both as a result of the merma and because of the contributions from the wood. Table 4 shows a high content of furfurals (5-hydroxymethylfurfural and furfural) in YPXFWs. These high values may be attributable to the Maillard reactions that occur during the heating of grape berries as they are directly exposed to sunlight to obtain the raisins. Moreover, a sugar medium with a pH of between 4 and 5 promotes these reactions [51]. As the wine was aged under an oxidative process in a Criaderas y Solera system, a significant increase in 5-hydroxymethylfurfural and furfural content was observed, again attributable to Maillard reactions, where the long aging period of these wines facilitates the contact between the highly concentrated sugars and the amino acids in the medium. 5-hydroxymethylfurfural is also yielded by the wood, to the thermal degradation of hexoses and rhamnose [14,52]. Furfural, on the other hand, has its origin in the wood itself, since it is produced from the dehydration and cyclization of pentoses during the thermal treatment the casks' wood is subjected to as part of its manufacturing process.

Evolution of Brown Color during Aging
The absorbance values measured were significantly different (ANOVA) depending on the aging scale. The absorbance at 470 nm allowed for the establishing of comparisons within the amber-mahogany-ebony range. The color of PXSWs increases with aging time (Table 4). This is due to hydroxylation, oxidation, and/or to the polymerization reactions of certain compounds in the wine that generally present a phenolic structure, such as quinones, among others, which change their appearance from colorless to colored [4,53]. The presence of certain metals, such as copper or iron, in wines has an influence on these oxidative processes, since they act as catalysts that facilitate the oxidation processes in PXSWs. Metal contents in PXSWs are higher than those determined in dry sherry wines of the same age [7,8]: 5-15 mg/L iron and less than 1 mg/L copper. All of the above explains the significant increments that have been registered when measuring the absorbance of these sweet sherry wines at 470 nm. Likewise, as the wines age, the transfer of certain compounds from the casks' wood also contributes to the wines' final color.

A Study on Predicting the Age of the Wines
A regression study has been proposed to determine those variables studied in the PXSW Criaderas y Solera system that presented a direct relationship with the age of the wine, such that a model could be developed to predict the age of the wines and to identify the most influential parameters that could facilitate wineries having a choice of routine analyses with which to control the quality of their wines.
As a first stage, a correlation study was carried out that included all of the variables in order to determine the degree of correlation between the predictive variables and age as the dependent variable. This study allowed for the removal potassium, isobutanol, and ethyl decanoate from the set of predictive variables, as they did not correlate significantly with the dependent variable, age, at a 95% confidence level. This first study also revealed that many of the remaining variables could be cross-correlated with each other, which would result in certain multicollinearity issues with respect to the final regression model.
In order to identify those variables with a greater weight in a potential regression model and verify whether the data have a latent structure that could be related to the age of the wines, a factor analysis using the remaining 46 variables was then conducted. A total of 84 instances of data for each parameter was evaluated. According to this analysis, 2 factors explained 95.531% of the variability, with VF1 accounting for 93.347% of the variance. In addition, according to the dispersion graph (VF1 vs. VF2), both factors are related to the age of the wine (Figure 2).
iance. In addition, according to the dispersion graph (VF1 vs. VF2), both factors are related to the age of the wine (Figure 2). Given this evidence, it was decided to continue the study with all of the variables whose VF1 coefficients had an absolute value equal to or greater than 0.5, which led to the removal of acetaldehyde, acetal, methanol, n-propanol, ethyl dodecanoate, calcium, nhexanol, and ferulic acid from the set of predictive variables. The next step consisted in searching for MLR models that would allow for the discrimination of those age-related variables that might be validated as markers or indicators of the wines' ages.
Given the large number of remaining predictive variables (39), a strategy that aimed at generating a model with the shortest possible number of independent variables was followed. For this purpose, a maximum number of five variables was established.
As a first step, the variables were grouped by families of compounds. The family variables include total acidity, volatile acidity, sulfate, total dry extract, sugar-free extract, absorbance at 470, and FCI. A "Multiple Regression" study was used to determine the best model for each group of variables. All of the predictive variables from zero to five were combined, and the statistics for their evaluation were then calculated, including the mean square error (MSE), the R 2 values adjusted to the degrees of freedom, Mallows' Cp, the p-values of the predictive variables and that of the model at 95% significance, and the Schwarz-Bayesian information criterion (SBIC). Thus, the best model was selected according to the following criteria: p-values < 0.05 of the predictor and the model variables; R 2 close to 100; and better (smaller) Mallows' Cp as well as SBIC criteria.
Once the best model for each group of variables had been obtained, the possibility of combining the different groups of variables was verified by following the same strategy. In this case, the only potential predictive variables to be considered were the variables that had been selected for the previous models.
In the case of the PX5 model, which combined all of the families of variables, and given that the procedure that had been applied for the previous ones had presented some Given this evidence, it was decided to continue the study with all of the variables whose VF1 coefficients had an absolute value equal to or greater than 0.5, which led to the removal of acetaldehyde, acetal, methanol, n-propanol, ethyl dodecanoate, calcium, n-hexanol, and ferulic acid from the set of predictive variables.
The next step consisted in searching for MLR models that would allow for the discrimination of those age-related variables that might be validated as markers or indicators of the wines' ages.
Given the large number of remaining predictive variables (39), a strategy that aimed at generating a model with the shortest possible number of independent variables was followed. For this purpose, a maximum number of five variables was established.
As a first step, the variables were grouped by families of compounds. The family variables include total acidity, volatile acidity, sulfate, total dry extract, sugar-free extract, absorbance at 470, and FCI. A "Multiple Regression" study was used to determine the best model for each group of variables. All of the predictive variables from zero to five were combined, and the statistics for their evaluation were then calculated, including the mean square error (MSE), the R 2 values adjusted to the degrees of freedom, Mallows' Cp, the p-values of the predictive variables and that of the model at 95% significance, and the Schwarz-Bayesian information criterion (SBIC). Thus, the best model was selected according to the following criteria: p-values < 0.05 of the predictor and the model variables; R 2 close to 100; and better (smaller) Mallows' Cp as well as SBIC criteria.
Once the best model for each group of variables had been obtained, the possibility of combining the different groups of variables was verified by following the same strategy. In this case, the only potential predictive variables to be considered were the variables that had been selected for the previous models.
In the case of the PX5 model, which combined all of the families of variables, and given that the procedure that had been applied for the previous ones had presented some difficulties, these issues were probably related to multicollinearity, as it generated models with appropriate selection criteria but where the p-values of the variables presented a lack of significance. Therefore, and in order to overcome this issue, this new model was generated using a backward stepwise selection procedure aiming to obtain an expression of the age of the wine based on eight variables. Table 5 presents the best models calculated for the PXSW samples. The PX5 model, which combines all of the previously selected variables, has a significance level above 95% and an R 2 (adjusted for GL) that explains 99.992% of the variance.  In order to verify the presence of multicollinearity and also to increase the stability of the model, a ridge regression was performed. This procedure forces the model coefficients towards zero, which reduces the risk of overfitting, the variance, any effects from the correlations between predictors, and the influence of the less relevant predictors on the model.
When this procedure was applied, the PX5 II model in Table 5 was obtained, which presents an R 2 adjusted to the degrees of freedom of 97.473 for a chain parameter of 0.065, showing variance inflation factors (VIF) coefficients below 3.
The resulting PX5 II model was influenced by the following compounds or analytical parameters: the increasing concentration with the aging of the wines of three of its most important phenolic acids, namely gallic acid, p-coumaric acid, and syringic acid, furfural, sugar-free extract, 2-phenylethanol, malic acid, and fertaric acid.
From the initial data matrix, six cases were extracted and used to validate the model after it had been built. Table 6 shows the results obtained.  Figure 3 contains the mean scores awarded by the panel to the samples evaluated, as well as the p-values resulting from the ANOVA applied to the three wines evaluated. Standard deviations of less than one for all of the wines and descriptors evaluated confirmed the high homogeneity of the panel. In general, nearly all of the descriptors (except for raisined fruit) demonstrated some discriminating capacity (p anova < 0.05) for the evaluated wines.

Tasting Sessions
Standard deviations of less than one for all of the wines and descriptors evaluated confirmed the high homogeneity of the panel. In general, nearly all of the descriptors (except for raisined fruit) demonstrated some discriminating capacity (panova < 0.05) for the evaluated wines.
The intensities corresponding to most of the olfactory and olfactory-gustatory descriptors increased after the aging of the YPXFWs, with the exception of raisined fruit, whose perception decreased somewhat (although not significantly) in comparison with the young wine. Nevertheless, it stayed above seven points in all cases. With regard to the sweetness of the wines, it was only with the over 30 years very long aging that its intensity increased to be statistically significant by more than 1 point and reached an average of 8.4 points. The rest of the PXSW descriptors increased proportionally with aging time. Thus, a clear increment in oak and toasted notes as well as a greater persistence and balance in the mouth were reported. Figure 3 allows a visual comparison of the sensory profiles of the wines evaluated. The multivariate processing of the data via factor analysis led to similar results. The factors were extracted based on the principal components. The factor analysis was conducted by applying a varimax rotation. The factors were selected according to an eigenvalue greater than one.
The analysis revealed 2 factors (F) that explained 90.4% of the total variability in the data (Figure 4a). F1 explains 77.3% of the variability; this factor is related to the aging time in the casks, since the ordering from left to right of the wines starts via the YPXFWs and The intensities corresponding to most of the olfactory and olfactory-gustatory descriptors increased after the aging of the YPXFWs, with the exception of raisined fruit, whose perception decreased somewhat (although not significantly) in comparison with the young wine. Nevertheless, it stayed above seven points in all cases. With regard to the sweetness of the wines, it was only with the over 30 years very long aging that its intensity increased to be statistically significant by more than 1 point and reached an average of 8.4 points. The rest of the PXSW descriptors increased proportionally with aging time. Thus, a clear increment in oak and toasted notes as well as a greater persistence and balance in the mouth were reported. Figure 3 allows a visual comparison of the sensory profiles of the wines evaluated.
The multivariate processing of the data via factor analysis led to similar results. The factors were extracted based on the principal components. The factor analysis was conducted by applying a varimax rotation. The factors were selected according to an eigenvalue greater than one.
The analysis revealed 2 factors (F) that explained 90.4% of the total variability in the data (Figure 4a). F1 explains 77.3% of the variability; this factor is related to the aging time in the casks, since the ordering from left to right of the wines starts via the YPXFWs and moves through the PX12YO and up to the PX+30YO. Figure 4b shows the weights corresponding to the sensory descriptors with regard to the first two factors of the space obtained, which should allow for the confirmation of their interpretation when discriminating between the profiles of the wines evaluated. Thus, a close correlation between F1 and the oak note can be observed, as expected, but one can also be observed between aromatic intensity, toast, and persistence and balance in the mouth, all of which have correlation coefficients with an F1 value greater than 0.9. These are the descriptors that increased with the aging of the wine. On the other hand, F2 explains 13.1% of the total variability in the data; this factor presents a strong negative correlation with the aroma of raisined fruit (r < −0.9), which explains the placement of the young PX wine at slightly lower positions than those of the aged wines. the oak note can be observed, as expected, but one can also be observed between aromatic intensity, toast, and persistence and balance in the mouth, all of which have correlation coefficients with an F1 value greater than 0.9. These are the descriptors that increased with the aging of the wine. On the other hand, F2 explains 13.1% of the total variability in the data; this factor presents a strong negative correlation with the aroma of raisined fruit (r < −0.9), which explains the placement of the young PX wine at slightly lower positions than those of the aged wines.

Conclusions
The characteristics of PXSWs, according to the analyses that have been conducted, are greatly influenced by the raisining processing of its grapes and by its oxidative aging. The decrease in alcohol content, due to the difficulty in maintaining the hermeticity of the casks, because of the high sugar content of the wine and the effect of the merma, causes many of the compounds studied to increase their concentrations during aging. Thus, an increase in the values of total acidity, volatile acidity, sulfates, glycerol, total dry extract, reducing substances, and sugar-free extract is observed as the PXSWs age. Tartaric acid content decreases with aging as it becomes insolubilized in the form of calcium and potassium salts, while the other organic acids analyzed-citric, malic, succinic, and lactic acid-increase their concentrations; of these, malic acid is the majority in the older wines. Among the volatile substances, it is worth highlighting the increase in the concentration of higher alcohols with aging, as well as esters: ethyl acetate, ethyl lactate, and diethyl succinate, which are also increased by the esterification carried out by their organic acids and the alcohol of the medium. On the contrary, a significant increase in fatty acid ethyl esters is not observed. Aging in oak barrels has shown a constant increase in FCI in all PXSWs, as well as in furanic derivatives as a consequence of Maillard reactions between sugars and nitrogenous substances in the wines. The combination of all of these factors leads to the remarkable organoleptic personality exhibited by this unique wine. With aging, the tonality in these wines changed from amber to ebony; additionally, the aromas and flavors, which in the YPXSWs were reminiscent of their origins (raisins, sweetness), changed to a greater complexity in the PXSWs that included the contribution of oak and toasted wood. This study has demonstrated the close correlation between the average age of PXSWs and the parameters that have been investigated through a multiple linear regression analysis. A model has been developed that has the capacity to determine the average age of these wines at over 97% confidence. A model of this kind should be considered as a powerful tool with which wineries can control wine quality, since it allows for

Conclusions
The characteristics of PXSWs, according to the analyses that have been conducted, are greatly influenced by the raisining processing of its grapes and by its oxidative aging. The decrease in alcohol content, due to the difficulty in maintaining the hermeticity of the casks, because of the high sugar content of the wine and the effect of the merma, causes many of the compounds studied to increase their concentrations during aging. Thus, an increase in the values of total acidity, volatile acidity, sulfates, glycerol, total dry extract, reducing substances, and sugar-free extract is observed as the PXSWs age. Tartaric acid content decreases with aging as it becomes insolubilized in the form of calcium and potassium salts, while the other organic acids analyzed-citric, malic, succinic, and lactic acid-increase their concentrations; of these, malic acid is the majority in the older wines. Among the volatile substances, it is worth highlighting the increase in the concentration of higher alcohols with aging, as well as esters: ethyl acetate, ethyl lactate, and diethyl succinate, which are also increased by the esterification carried out by their organic acids and the alcohol of the medium. On the contrary, a significant increase in fatty acid ethyl esters is not observed. Aging in oak barrels has shown a constant increase in FCI in all PXSWs, as well as in furanic derivatives as a consequence of Maillard reactions between sugars and nitrogenous substances in the wines. The combination of all of these factors leads to the remarkable organoleptic personality exhibited by this unique wine. With aging, the tonality in these wines changed from amber to ebony; additionally, the aromas and flavors, which in the YPXSWs were reminiscent of their origins (raisins, sweetness), changed to a greater complexity in the PXSWs that included the contribution of oak and toasted wood. This study has demonstrated the close correlation between the average age of PXSWs and the parameters that have been investigated through a multiple linear regression analysis. A model has been developed that has the capacity to determine the average age of these wines at over 97% confidence. A model of this kind should be considered as a powerful tool with which wineries can control wine quality, since it allows for the determination of a wine's age based on just eight parameters: syringic acid, p-coumaric acid, sugar-free extract, malic acid, 2-phenylethanol, gallic acid, fertaric acid, and furfural.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.